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Restricting the functional integral to the Gribov region Q, leads to a deep modification of the behavior 
of Euclidean Yang-Mills theories in the infrared region. For example, a gluon propagator of the 

,2 

Gribov type, ^4^-_4 , can be viewed as a propagating pair of unphysical modes, called here i-particles, 

with complex masses ±47^. From this viewpoint, gluons are unphysical and one can see them as 
being confined. We introduce a simple toy model describing how a suitable set of composite operators 
can be constructed out of i-particles whose correlation functions exhibit only real branch cuts, with 
associated positive spectral density. These composite operators can thus be called physical and are 
the toy analogy of glueballs in the Gribov-Zwanziger theory. 

PACS numbers: ll.lO.-z, 03.70.+k, 11.15.-q 



I. INTRODUCTION 



A very peculiar feature of the strong interaction is the fact that quarks and gluons cannot be observed as free par- 
ticles, i.e. they are absent from the physical spectrum. As the theory of quantum chromodynamics (QCD) describes 
the strong interaction, it should account for this property. A first step towards a full picture is the investigation of 
the gluonic sector alone, thus ignoring contributions from quarks. Even this simpler case of Yang-Mills theory is very 
intricate and is widely studied from both analytical and numerical points of view. Its spectrum is believed to consist 
of so-called glueballs, color singlets made up by gluons only^. Yang-Mills theory then has to explain why gluons do 
not exist as asymptotic states, whereas glueballs do. 

A possible and well investigated mechanism for gluon confinement is provided by the Gribov-Zwanziger framework 
[21 [3], which enables us to take into account the effect of Gribov copies [2]. Although the existence of Gribov copies 
is a general feature of the gauge fixing procedure [4,, we shall refer here to the case of the Landau gauge, dfj^A^^ = 0, 
which is a renormalizable gauge in the continuum, while also possessing a lattice formulation. The problem of gauge 
copies is remedied by restricting the domain of integration in the Feynman path integral to the Gribov region Q, 
defined as the set of field configurations which obey the Landau condition and for which the Faddeev-Popov operator, 
M"'' = -d^Df = -{S'^^d'^ + .g/"^M^9^), is strictly positive, namely 

n^{A; d^A"^ = 0, X"'' > } . (1) 

The region fl is convex and bounded in all directions in field space, and every gauge orbit passes through fl [5J. 
Its boundary d^l is known as the first Gribov horizon, where the lowest eigenvalue of the Faddeev-Popov operator 
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vanishes. As shown in [3", 'S', the restriction to the region n is achieved by adding to the Yang-Mills action the so-called 
horizon term 

Sh =7'/ d'xh{x) =7'/ d\ J d^yDl^{x){M-'r\x,v)Dl^{v) , (2) 

with 

Df = + (3) 

the covariant derivative. The action thus becomes 

S = Sym + Sgf + Sfi , (4) 
where the Yang-Mills and gauge fixing part of S are given by 

Sym = 7 / d^xF^^F";, , (5) 
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The field 6" is the Lagrange multiplier enforcing the Landau gauge condition, and {c°',c°') are the Faddeev-Popov 
ghosts. The massive parameter 7, known as the Gribov parameter, is not free, being determined in a self-consistent 
way by a gap equation called the horizon condition f3^, namely 

{hix))^diN'-l), (7) 

where N is the number of colors and d the space-time dimension. 

Although being nonlocal, the horizon term can be cast in local form by introducing a suitable set of auxil- 
iary fields {'f'j^, (f'^'', w"'', w^''). The fields {if'^, (p^f') are a pair of complex conjugate bosonic fields, whereas (oj"'', w"^) 
are anticommuting. Notice that these fields carry two color indices and one Lorentz index. This is necessary in order 
to account for the required number of degrees of freedom. Thus, the localized form of the action reads 

Sgz = Sym + Sgf + Sioc , (8) 

where 

Sioc^ J d^x{^l^M^\l^-Cjl^M'^''^l^-gr''^{d,Cjf){Dl^c^)^^^^ . 

(9) 

The last term is introduced in order to be able to rewrite the horizon condition ([t]) as 

5T , , 

^=0, (10) 

where F is the vacuum energy, i.e. 

= / e--^-^^ , (11) 



and [d^] stands for the functional integration over all fields. In terms of the auxiliary fields {ip1^,'p'^), the horizon 
condition (10 1 takes the form 

(sr^^A^^ - ^^^)) = 2d{N^ - 1)7' . (12) 

Remarkably, the action ([8| turns out to be multiplicatively renormalizable to all orders, thanks to the existence of a 
rich set of Ward identities which can be established at the quantum level [51171-fT^. Let us mention in particular that, 
unlike the case of the Faddeev-Popov action, the Gribov-Zwanziger action does not enjoy exact BRST symmetry, 
which turns out to be softly broken by terms proportional to the Gribov parameter 7 |12f[I^ . Nevertheless, the 
breaking term can be kept under control at the quantum level, giving rise to softly broken Slavnov-Taylor identities 
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which constrain very much the form of the possible counterterm. We also refer to fT5' for a few peculiarities 
concerning the renormalization analysis of the Gribov-Zwanziger theory. 

The restriction to the Gribov region Q entails deep modifications in the behavior of the theory in the infrared region. 
Let us give a look, for example, at the gluon propagator stemming from the action (|8| 

{A'^,ik)Ali-k)) S^" (j,. - ^) ^ , 7^ = AT/ . (13) 



Expression (13 1 is suppressed in the infrared region, while displaying complex poles at k — ii-y . This structure 
does not allow us to attach the usual particle meaning to the gluon propagator, invalidating the interpretation 
of gluons as excitations of the physical spectrum. In other words, gluons cannot be considered as part of the 
physical spectrum. In this sense, they are confined by the Gribov horizon, whose presence is encoded in the explicit 



dependence of expression ( 13 1 on the Gribov parameter 7. Another way to see that gluons are not physical particles is 



the observation that the propagator (13 1 has negative norm contributions and thus it violates positivity |16j. Since a 



positive-definite norm is required for a probabilistic interpretation |17|[T5] . we conclude that the excitations described 



by eq. ( 13 1 are confined. One might argue that further nonperturbative effects besides the Gribov horizon might 



conspire to change this picture. However, qualitatively one can obtain the same result with functional methods |,19,-.24] . 

Even if the restriction to the Gribov region, with the consequence of an infrared soft breaking of the BRST 
symmetry, provides a natural mechanism for gluon confinement, one is left with the nontrivial task of constructing 
a set of physical composite operators made of the fields entering the Gribov-Zwanziger action ([s]), whose correlation 
functions can be given a physical meaning in Minkowskian space. A possible strategy to work out this kind of analysis 
is to look at the spectral properties of the correlation functions. In practice, this amounts to face the challenging 
task of constructing a suitable set of composite operators whose correlation functions can be cast in the form of a 
Kallen-Lehmann representation, i.e. a spectral representation with a positive spectral function, and whose analytic 
continuation in the complex Euclidean /c^-plane exhibits a cut along the negative real axis only. Such a spectral 
representation would thus imply that, when moving to Minkowskian space, the cuts are located along the positive 
real axis. Moreover, positivity of the spectral function then guarantees that a meaningful interpretation in terms of 
states of a physical spectrum can be attached to those operators. This is precisely what one would expect from a 



confining theory. This is a highly nontrivial task, given the complexity of the Gribov type propagator ( 13 1 as well as 
of the Gribov-Zwanziger action ([8|. 

The first step towards a detailed study of the analyticity properties of the correlation functions within the 
Gribov-Zwanziger framework was undertaken by one of the authors [3], who evaluated at one loop order the 
correlation function of two gauge invariant operators, namely 

where F'^{x) — F'jl^j{x)F^^^{x). The results found in [3^ can be summarized as follows: 

G'(fc2) = G^'^^^k^) + G""P'^y"(fc2) , (15) 

The unphysical part, G^'^'^^^^ {k"^) , displays cuts along the imaginary axes beginning at the unphysical values 
k^ = zLAij"^, whereas the physical part, G^^y^{k^), has a cut beginning at the physical threshold k^ = —27^. 
Moreover, the spectral function of G'p^^''(A;^) turns out to be positive, so that it possesses a Kallen-Lehmann 
representation [3]. As such, G^^^^{k^) is an acceptable correlation function for physical glueball excitations. What 
is interesting in expression ( |15[ ) is that a physical cut has emerged in the correlation function of a gauge invariant 
quantity, even if it has been evaluated with a gluon propagator exhibiting only unphysical complex poles. 

Obviously, the big challenge lying ahead of us is to find out if one can get rid of the unwanted unphysical 
cut. We should like to mention an alternative possibility which at the moment remains at the level of a scenario, 
in fact, a scenario with a deus ex machina: it is possible that physically satisfactory analyticity properties emerge 



only when the horizon condition, eq.(10l, is satisfied. This is difficult to investigate because the horizon condition is 
nonperturbative, and we can offer no evidence in support of this happy eventuality. In the present investigation we 
are limited to perturbative calculations. 



The present paper aims at pursuing such an analysis. To go further, we should figure out a plausible mecha- 
nism which deforms local gauge invariant operators Oi such that their expectation values have only physical cuts. 
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This deformation can be done by adding to Oi local quantities depending on the fields of the Gribov-Zwanziger 
action. As it will be shown, this deforming mechanism will allow, in particular, to get rid of the unphysical cuts in 



the lowest non trivial order for the above-mentioned two point function, eq.(14l. The possibility that the mechanism 
can be generalized for other operators, and at any given order of perturbation theory, could be envisaged as follows. 
One first considers the usual observables of Yang-Mills theory: basically the elements of the BRST cohomology, i.e. 
the color singlet local gauge invariant operators built up of the field strength F^^, and its covariant derivatives. For 
a given canonical dimension, one has generally a mixing that is induced by renormalization |25l 126]. For instance, 
an operator such as F'^ = f^^F^^^FlpF^^ can mix with other BRST nontrivial operators such as F^^{D'^F^^)°- , etc. 
Moreover, the softly broken Slavnov- Taylor identities indicate that, besides the expected mixing with BRST exact 
operators, also mixings with operators which are not BRST invariant will, in general, occur. One hopes, however, 
that only suitable combinations of all these operators can be defined, such that their expectation values only exhibit 
physical cuts, with mutual compensation of all terms with unphysical cuts. One is thus led to a kind of boot-strap 
mechanism, where one starts with a particular mixing at the tree level, such that the absence of unphysical cuts be 
enforced order by order in perturbation theory. In other words, at the tree level, an observable should be written in 
the following form 

Q ^ ^ ^ Q:^((^'^)C'^^^^^ invariant) ^ (^g'^^C)^'^°'^ BRST invariant) ^ (16) 

i 

The values of the coefficients a*(0) and /3*(0) at their lowest order are supposed to be determined by the requirement 
of having only physical cuts. Our point of view is thus as follows. The BRST symmetry characterizes the 
action and the observables of the theory in the perturbative short distance regime as elements of its cohomology. 
This symmetry is however softly broken in the long distance regime. Even if the presence of the soft breaking 
enables us to determine a renormalizable local action under the form of the Gribov-Zwanziger action, it turns 



out that the local observables have to be modified by the addition of suitable terms as in eq.(16l, where the 
positivity of the corresponding correlation functions is required to hold. When all the auxiliary fields are integrated 
out, this is tantamount to modifying the observables and the action by nonlocal terms, proportional to the 
nonperturbative Gribov parameter 7. The remaining of this paper will consist in providing a certain number of ex- 
amples in order to check the plausibility of this framework, by using the predictive power of local quantum field theory. 

Our main purpose here is therefore that of investigating how the use of a confining Gribov type propagator 
can allow us to obtain examples of correlation functions exhibiting real cuts only and positive spectral functions. Our 
analysis relies on the introduction of what we shall call j-particles: a pair of fields with complex conjugate masses 
which emerge in a natural way when dealing with a Gribov type propagator. In fact 



1 1/1 1 



(17) 



from which it is apparent that a Gribov propagator can be associated to the propagation of excitations with complex 
conjugate masses ±17^. As it will be discussed in details, the use of z-particles enables us to provide examples of 
operators whose correlation functions exhibit the desired analyticity properties. More specifically, we shall be able to 
show that local composite operators built up with pairs of i-particles display cuts along the negative real axis in the 
complex Euclidean fc^-plane, while giving rise to positive spectral functions. 

The paper is organized as follows. In Sections [ll]|VI| we discuss our ideas with a toy model, which is intro- 
duced in Section |TT] using a scalar field possessing a confining propagator of the Gribov type. This will enable us to 
illustrate in a simple way how the action of the toy model can be cast in diagonal form through the introduction 
of i-particles. Section |III| is devoted to a detailed study of a first example of a local operator, constructed from 
i-particles, whose one loop correlation function possesses a Kallen-Lehmann spectral representation. Section |IV| 
provides further details of the calculation of this correlation function in the complex fc^-plane. In Section |v] we 
face the evaluation of correlation functions with several loops. The evaluation of the spectral representation for 
two and three loop integrals obtained from operators built up with i-particles will be worked out. In particular, it 
will be shown that the spectral function at higher order can be obtained by a kind of convolution of the spectral 
functions obtained at lower orders. This provides a useful iterative procedure, allowing us to generalize the argument 
to the n-loop case. In Section |VI| we discuss the meaning of the i-fields within the toy model. We argue that 
those fields might be regarded as rather useful variables in order to construct sensible composite operators within a 
local quantum field theory framework. In Section |VII| we return to the Gribov-Zwanziger action. We discuss how 
i-particles emerge and how they can be employed to give examples of one loop correlation functions which, unlike 



expression (15 1, display cuts along the negative real axis only while having positive spectral functions. Even if a 
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systematic construction of a set of meaningful operators in the presence of the Gribov horizon is far form being 
realized, the introduction of the i-particles might be seen as an interesting path in order to provide at least some 
examples of operators displaying good analyticity properties, something which can already be regarded as a nontrivial 
achievement. In the Conclusions a few points for future investigations will be spotted. Finally, several useful details 
about the evaluation of the spectral densities are reported in the Appendices. 



II. A SCALAR FIELD THEORY TOY MODEL 



A. Constructing the toy model 

A simple way of constructing a field theory model exhibiting a confining Gribov type propagator is through a scalar 
field -0 whose Euclidean action is nonlocal, being specified by 

S^J d^x^-^j(^~d' + 2^^^. (18) 

The massive parameter 9 is introduced by hand and is akin to the Gribov parameter 7. As it is easily seen, the 
resulting propagator is in fact of the Gribov type 

WkM'k)) = . (19) 



Although being nonlocal, expression ( 18 1 can be cast in local form by introducing a pair of bosonic complex conjugate 
fields {<p,(p) and a pair of anticommuting fields (tl),a;), so that 

5 = J d-^x (^^tP{-d^)i: + <f{-d^)^ + e^ij{^-'^)-u!{-d^)u?j . (20) 

Let us have a look at the propagators of the model. To that purpose, we introduce new variables (J7, V) in order to 
achieve a partial diagonalization of the action 

U + iV 

In terms of the fields ([/, V^), the action of the toy model reads 



(22) 



From the above expression, the propagators can be worked out 

{VV) 







+ 26*4 








p' 




+ 26*4 


— ? 






+ 26*4 


1 




p2 





(UU) = ^ . (23) 

Having evaluated the propagators of the fields {ip, U, V), we can now check what the propagators in terms of the fields 
{f, If) are. One finds 
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+ 26*4 



(24) 



B. Introducing the i-particles 



Expression ( 22 1 can be cast in complete diagonal form by making a further change of variables: 



so that 



(25) 
(26) 

(27) 



From this expression one immediately sees that the fields A and rj correspond to the propagation of unphysical modes 
with complex masses ±i\/26l2. These are the i-particles of the model, namely 



(A(fc)A(-fc)) 
{ri{k)v{-k)) 



It is important to notice here that requiring 



At 



1 

1 

fc2 - i\/26l2 



7^ = -A 



(28) 
(29) 

(30) 



will ensure that (27 1 corresponds to a Hermitian action. 



III. CONSTRUCTION OF COMPOSITE OPERATORS AT ONE LOOP 



A. Preliminaries 



Having cast the action in diagonal form, we face now the construction of a set of composite operators made 
out of z-particles, whose correlation functions exhibit real cuts only. As the Gribov propagator ( 19 1 and the di- 
agonal form of the action ([6| suggest, we shall look at composite operators constructed by means of pairs of i-particles. 

The simplest example which one can consider at one loop is that of the dimension two composite operator 
consisting of one i-particle of the type A and one i-particle of the type r/, namely 



(31) 



which is indeed a real (Hermitian) operator, according to the prescription (30) 



The correlation function {Oi{k)Oi{—k)) in d Euclidean dimensions is given by 

d'^p 1 1 



(Oi(fc)Oi(-fc)) 



{2tiY [k - pY - iV^e^ p2 l^Q2 



(32) 



^ There are more propagators than the ones shown, but these are not relevant for the calculation presented. 
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By direct inspection of the action (|6| , it follows that the correlation function of three operators Oi (x) vanishes 

(Oi(x)Oi(y)Oi(z)> =0. (33) 
Only correlation functions with an even number of operators Oi are nonvanishing. 

The complex nature of the i-particle masses may raise concerns about the consistency of the theory. A press- 
ing question is that of the unitarity. One has to be sure that these correlation functions have a well defined 
probabilistic interpretation. In order to settle this question, in what follows we shall provide an explicit evaluation of 
the spectral function associated with the Kallen-Lehmann representation of the two-point function (Oi(fc)Oi(— fc)). 

B. Kallen-Lehmann representation in d = 4 



We start our analysis by evaluating the integral ( 32 1 . Our way to proceed is that of casting it into the so-called 



spectral representation, which has a precise and powerful meaning, both in complex analysis and quantum field theory. 
It is given by 



(Oi(fc)Oi(-fc)) = / drpir) 



1 



fc2 



(34) 



where the quantity tq > stands for the threshold. Once expression (34) has been obtained, we take it as our very 
starting point. We introduce the complex function 



(35) 



Therefore, from complex analysis, it follows that F{z) is an analytic function in the cut complex plane, where 
the interval (— oo, — tq) has been excluded. Moreover, when moving from Euclidean to Minkowskian space, i.e. 



^Eucl 



^Mink' expression (35) gives the spectral representation of a quantity exhibiting a discontinuity along the 



positive real axis, starting at the threshold tq and extending till +oo. One requires also that the function /o(t) is 
positive in order to have a probabilistic interpretation. Indeed, using Cauchy's theorem, it is a relatively easy exercise 
to show that 



p(t) = lim \F{—T — ie) — F(~t + ie)] , with t > tq . 

2m e-i-o+ 



(36) 



i.e. p{t) is directly related to the discontinuity of the propagator F{z) along its branch cut. Via the optical theorem, 
we then also know that p(t) is proportional to a cross section, which evidently needs to be positive in a meaningful 
quantum field theory. 

In order to obtain the spectral representation, we employ Feynman parameters in the general d-dimensional 
situation in order to combine the denominators in ( 32 1 



{Oi(k)0^{~k)) 



1 



1 



d'^p , 1 



(2^)' JO 



dx- 



dx- 



[a;(fc2 -^2p-k- 2iV2e^) + p^ + iV2e^] ' 
1 



[g2 + {x- a;2)fc2 - {2x - l)iV20^ 
where we have defined q = p — kx. Using now the identity 

d^q 1 _ 1 r(n-f)^^2^._„^ 



{2nr [g2 + A2]" (4^)f r(n) 
with n = 2 and A'^ = {x - x^)P - {2x - l)iV29'^, we obtain 



(Oi(fc)Oi(-A:)) = 



(4^) 



dx 



{x - x^)k^ - {2x - l)iV2e^' 



#-2 



(37) 



(38) 



(39) 
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We now consider the case of d = 4. We define 



= (Oi(fc)Oi(-fc)) 



(40) 



We start from (39l and act on it with gla- This regularizes the original integral for F(fc^), which is ultraviolet 
divergent. After setting d = 4, we find 



dx- 



x{\ — x) 



■{l-x)B +xi-i/2 



167r2 



dx 



x{l — x) 



x{l — x)k^ — xi + i/2 ' 



(41) 



where we temporarily switched to units 2\/29^ = 1 for notational convenience. The substitution s 



^ _ \+s+^^\+ ^ ^ brings us to 



.F(fc2) 



1 



167r2 
1 

167r2 



+C30 



ds- 



oo 

+ 00 



ds 



(fc2 



ids -1 + s + \/l + 



(42) 



where we employed partial integration. There is no problem at s = if fc^ > 0, since lims_>.o ^+s^^\ ±s^ _ xhere 
are no poles in the upper half s-plane for fc^ > 0, so we can deform the contour to be located around the cut for 
s S [zcxDji]. Setting s — ir, we compute (42) as 



d 



F(fc2) 



1 



167r2 
1 



~dT 



IT — l\J T 



(fc2 +t)2 



2iT 



167r2 



-dr 



-1 + ir + r2 



(fc2 +t)2 



2iT 



167r2 



(F 



We can subsequently integrate this expression from to A:2^ finding 



F(fc2)_F(0) 



1 



167r2 



1 



+ fc2 



- Mt. 



or, by restoring the units, 



F{k'^)~F{Q) 



167r2 



8614 



fc2 



-Mr. 



(43) 



(44) 



(45) 



From this expression, we notice the importance of the subtraction of F(Q) to find a finite result, otherwise we would 
find a divergent spectral integral. The spectral density can be read off from ( 45 1 , 



P{r) 



1 Vt2 - 86*4 



167r2 



(46) 



Interpreting this result as an expression of the physical spectrum of the theory, we see that it has a threshold at 
2\/26i2. The same calculation can be repeated in the case of two particles with real mass /i, where the spectrum is 
found to begin at the threshold 4/i2 = (y^ -f yi)2 the present case we have an analogous situation for the i-particles. 



Note also that the spectral function is 



even though they have complex masses: 2^/29^ ~ ^e'4 24( 

positive in the range of integration. This is a very important feature for the theory to have a physically meaningful 
probabilistic interpretation. 



Finally, an explicit integration of ( 44 1 leads to 



F{k^)-F{{)) 



1 



167r2 



TT 
2fc2 



fc2 



arccos(fc^ 



(47) 



In FIG. 1(a) and 1(b) we have displayed the (rescaled) real and imaginary part of F{k^). The cut for z G [— oo, — 1] 
is clearly visible. In Appendix |A] we have determined the spectral density p(r) for the case that d = 2 or d = 3 while 
verifying the d — 4 result, based on two quite different calculational schemes. Also in the cases d = 2,3, a positive 
spectral density is found. 
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G,(t-) 




FIG. 1: Plots of Gfl(fc2) = 167r2Re[F(fc2) - F(0)] and Gi{k'^) = 16nHm[F{k'^) - F(0)], respectively, with F{k^) - F{0) given 
in ||47|. 



IV. A CLOSER LOOK AT THE ANALYTIC CONTINUATION BY MEANS OF THE SPECTRAL 

REPRESENTATION 

A. The case of complex masses 



As already explained in Section III the knowledge of the spectral density p{t) allows us to obtain the analytic 
continuation of the correlation function for complex k^, whereby it not necessarily holds that fc^ > 0. The main 



formula is given by ( 35 1 , by means of which we are able to define the analytic continuation of the original momentum 
space integral (32 1, which gives for the correlation function'^, 



1 

4^ 



1 



(48) 



For any real external momentum vector fc, so that fc^ > 0, the integral (48 1 certainly makes sense. For such vectors, 
we can reexpress ( 48 ) as 



1 



1 



dx 



47r2 7 ^F^-V2p2 + p_2p-fc + i/2 ^'^ h xk^ ~ x^k^ + ix - i/2 ' 



(49) 



as the Feynman (or Schwinger) parametrization trick (39 1 is valid for fc^ > 0. We have again set 2y/29^ = 1 



The momentum integral (48) is however also well-defined for a more general set of k. More precisely, even 
for complex vectors k, (48 1 can be still well-defined, without having that F{k'^) given by (49) is a valid representation 



of it, since it is possible that neither the Feynman nor Schwinger trick applies for such vectors k. In Appendix |a] we 
have calculated the equivalent of (46) for d = 2, see eq. (A8l. Just as in the previous section, the spectral density 
Pd=2{T) was obtained using fc^ > 0. It is therefore nontrivial to verify that 0{k^) — Fi{k'^) for any choice of k 
whereby 0{k^) exists, with 



dr 



1 arccos(fc2) 



(50) 



which is obtained using the spectral representation (A8|. 



^ We shall work with d = 2 here for simplicity. 
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That this constitutes a not so trivial check, becomes more apparent when one reaUzes that once having ob- 
tained (49 1 for fc^ > 0, one can also use the r.h.s. of ( [49| to define a function for any complex value of fc^. Indeed, 
the integration can be done exactly, by making use of |27j 



dx 



- In " ^ ^ 
b a b 



(51) 



valid for any complex number a and b. The ill-definedness of the integral in the l.h.s. of (51 1 for G [0, 1] corresponds 
exactly to the branch cut of the In in the r.h.s. of (51 1. We can rewrite the integral (49 1 in the following form 



1 

An 



dx 



xk^ 



x'^k^ + ix ~ i/2 



1 

47r 



dx 



1 



Vk^ - 1 



1 



1 



by introducing partial fractions, with x± = ^'^^ ^t"^^ — ^- Subsequently, using (51 1 twice, yields 



F^(k^ 



47r - 1 



In 



In 



(52) 



(53) 



For fc^ > 0, one verifies that Fi{k^) = F2(fc^), but the situation changes drastically in the complex fc^-plane. For 
Re(fc^) > 0, it stills holds that -F'i(fc^) = ^2(^^)1 but the functions differ for Re(fc^) < 0. This difference can be 
traced back to the branch cut of F2{k'^), which is given by the complete imaginary axis. This observation follows 
immediately as the integral (52 1 is well-defined for any fc^ G C, except for k^ = is, s E R. In sharp contrast, Fi(fc^) 
displays a cut located on the negative real axis, as it was already explained before. 

We shall now argue that only Fi(fc^) obtained via the spectral representation gives a decent analytic continu- 
ation of the original momentum integral 0{k^). Let us thus choose an external momentum vector k which can be 
complex, or 

k = Re(fc) -I- ilm(fc) = kii + ikj . (54) 
We shall restrict ourselves to complex momenta of the type k = {u + iv, 0), which allow to reach any value of 



2iuv . 



(55) 



Although we already chose fc^ to be the argument of 0{k^), this is not a priori clear when k is not real. However, we 



notice that the integral in the r.h.s. of (48), if it at least exists , will define an analytic function of k 



■ iv, which 



can be easily checked by means of the Cauchy-Riemann equations. Since we still have the 0(2) rotational symmetry, 
as we can rotate fc^ and kj simultaneously, which thereby defines a real rotation of the complex vector k — kcx, it 
appears that ( 48 ) must be a function of fc^ . 



We then rewrite ( 48 1 as 



1 

4^ 



dpxdp. 



^ Px + Py — i/2 p"^ + p'l + u'^ — v"^ + 2iuv — 2px{u -\- iv) + i/2 



(56) 



One could wonder whether it would be possible to obtain finer results if we extend our choice of fc? Namely, using 
k — {u + iv,u' + iv'), it is not unconceivable that one could do better. To refute this possibility, we notice that we 
can already choose without loss of generality k = (u + iv,u'), since we can always rotate our basis to bring k into this 
form. The corresponding integral would read 



0'(fc2) 



1 



dpxdp. 



1 



^ Px + Py — i/2 p"^ + p'^ + — v"^ + 2iuv + — 2px{u + iv) — 2pyu' + i/2^ 



(57) 



but the simple shift py — 
without loss of generality. 



Py + u gives 0{k^) — 0'{k^), from which we conclude that we can set k — {u + iv,0) 



* We shall shortly see that such k do exist. 
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To proceed, the integral (56) is well-defined if the integrand is free of poles for p e 
when such poles might appear. The only possibility occurs when^ 



Px 



1 

Av 



Let us thus check 



(58) 



which is necessary to kill the imaginary part in the second denominator, and consequently if 

1 



(59) 



Hence, if we take \v\ < 1/2, with u arbitrary, the integrand of (56) is a regular function of and py, and the integral 
will exist as such. It is interesting to notice that this does not mean that the Schwinger and/or Feynman trick is 
applicable for any such vector k. Nevertheless, the integral has a well-defined value. As such, it must coincide with 
other ways to compute or define it. 

We can now discriminate between Fi(fc^) and F2(fc^) by taking a test value in the left half plane, such that 
the original integral (56) exists. Its value should then coincide with either Fi{k'^) or F2{k'^). Actually, we can 



immediately motivate that F2{k ) cannot be correct in the whole complex plane. Looking at ( |55[ ), it is clear that 
/c^ e iM if w = ±w, in which case we have fc^ = ±2iv'^. But as we have shown, for \v\ < 1/2, the original integral 
is well-defined, meaning that there cannot be a cut on the whole imaginary axis, whereas F2{k'^) has a cut for any 
k^ £ iR. This shows that F2(fc^) cannot be the analytic continuation of the original momentum integral. In fact, 
Fi{z) also describes the analytical continuation of F2(z) with Re(z) > to the left half complex plane, defined 
by Im(z) < 0. Since the whole imaginary axis is a cut of F2{z), it is impossible to find an open region of C 
where -F2(-z)im(z)>o and -F2(-Z2)im(z)<0) coincide. Said otherwise, -F2(^)im(z)>o and -F2(^2)im(z)<o are not analytically 
connected. 



As a second illustration, let us look at the test value 

(u*, w*) 

It can be checked that in this case, we have 

Fi{kl) « 0.254 - 0.028i, F2ik^) « -0.239 - 0.024i . 



1 1 

4' 3 



-7 i 
144 6 



We can also calculate ( 56 1 directly. We write 



Oiki) 



47r2 



dpy 



+ 00 



dpx 



pI +pI - if^Px +Py + ul — vl + 2m*i;* — 2px{u^, -|- w*) + i/2 



(60) 



(61) 



(62) 



and we use the residue theorem for the paj-integration. We close the contour along the upper hemisphere at cxd, where 
the function vanishes, and we notice that 2 of the 4 poles are located within the contour, being 



Pxl = 



i-2pl 



V2 



Px2 = 



\:(3 + 4i- 6V2 J-i - 2pl) , 
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(63) 



which is independent of py. So, 



O(fc^) = 
fiPy) - 



36%/2 



+ CX3 



dPyfiPy) 



i - 2pl [7 + 120i + (36 + 48i) J-2i - Apl M - 2pl [7 - 168^ + (36 + 48^) J2i - Apl 



We can exclude here the v = case as this corresponds to the anyhow well-defined case of real external momentum. 
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(64) 

Since py G M, the remaining integral is well-defined^ and can be computed using direct primitivation. After the smoke 
clears, one finds ^ 



0(fc2) w 0.254 - 0.028i, 
a value which nicely coincides with -F'i(fcJ), while differing from F2{k^). 

We have provided a nontrivial verification of the correctness of the analytically continued function i^i(fc^) 
fact, much more can be learned from our analysis. Let us take expression (551 into reconsideration. 



2iuv . 



(65) 



In 



(66) 



As we have shown, problems can only occur for > 1/2. But even restricting ourselves to < 1/2, eq. (66 1 means 
that we can still reach a relatively large region of the whole complex /c^-plane without encountering difficulties, since 
u is arbitrary. For example, for m = 0, we can end up on the negative real axis, but there can be a cut on the negative 
real axis only for fc^ < — 1/4. As another example, we repeat that by taking u = ±v, we conclude that a cut on the 
imaginary axis is a priori possible, but only when |Im(A:^)| > 1/2. 



In FIG. 2, the shaded parabolic region corresponds to the values of fc e C for which the integral (56 1 is cer- 
tainly well-defined. We clearly observe a small region in the left half plane where F{k^) should be well-defined, even 
without invoking analytical continuation. The results obtained do, of course, not mean that there must be a cut for 




FIG. 2: Parabolic region in the complex plane where the original momentum integral 0{k^) exists. 

fc^ S [— c», — 1/4] or that there is a "cut region" outside of the displayed parabola. But we have discovered a certain 
region where the analytic continuation has to be consistent with the value of the integral. Certainly, the analytic 
continuation of Fi(fc^) is consistent with the value of the integral within the boundaries of FIG. 2. 

To close this section, we mention that we could also have started with the alternative definition 



Oik' 



47r2 



d2-HP_2p'.^_j/2F^+i/2 ' 



(67) 



® We will never come into the vicinity of the cut of the occurring roots. 
^ We did not write the closed expression, as it is rather intransparent. 
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which corresponds to switching the momenta running in the two legs. For real external momenta k, we obviously 
have 0(k^) — 0{k^), but for complex fc, we can no longer perform a translation on the real integration momentum p 
to prove this. 



However, again setting fc = (u + iv, 0), it is easily checked that 
is well-defined, while from FIG. |3(a)[ |3(b)[ \A{a._ " 



< 1/2 is sufficient to also guarantee that 
4(b} it is clear that 0{k^) = 0{k^) for u arbitrary, \v\ < 1/2. 
For completeness, we have also shown i^i(fc^) in FIG. 3(c) 4(c) which illustrates that indeed O(fc^) = O(fc^) = Fi{k'^) 
over the parabolic /c^-region shown in FIG. 2. To avoid confusion, we point out that we have plotted 0(u, v) = 0{k^), 
with fc^ = — f ^ + 2iuv, and analogously for the other functions. We notice that this is equivalent with having the 
Schwartz refiection principle, 0(fc^''') = 0{k'^)\ a property shared with Fi(fc^). 






FIG. 3: Re(0), Re(0) and Re(fi) 






FIG. 4: lm(0), lm(0) and Im(Fi). 
Although we focused exclusively on the d — 2 case in this section, similar conclusions can be drawn for c? = 3 or c? = 4. 



B. The case of two real masses 



In order to corroborate the previous nontrivial conclusions about a relatively simple Feynman integral with complex 
masses, we find it instructive to also include a similar analysis of the probably more familiar case of two real and 
positive masses, being mi and m2- This has been investigated in great detail in e.g. [5S], albeit in Minkowskian space. 



The analog of (48 1 is given by the following correlation function. 



02{k') = A / d^P ^ ^ 9 = ^ = ■ (68) 

Translating the results of [55] to Euclidean space, it was shown that this function of k^ has a positive spectral density 
in combination with a branch cut on the negative real fc^-axis starting from —(mi + TO2)^ until — oo. 
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Let us now also investigate the region where this integral actually exists. We can again restrict ourselves to 



complex momenta of the type k ^ {u + iv, 0) without loss of generality. Inserting this in eq. (68 1, we obtain, 



_±_ f^p^^p 1 1 

47^2 J ^ y p'2 _^ p2 j^2 p2 ^ p2 _|_ y2 _ ^2 _|_ 2iuv — 2px {u + iv) + ni^ 



Let us check when poles can emerge. The first denominator is always positive, however, the second denominator can 
have poles when the imaginary part vanishes. 



2iuv — 2ipxV = 0<=>f = or u ~ px , 
simultaneously with the real part being equal to zero, 



Inserting the first solution of eq. (70), i.e. w = 0, in the equation above, we find 



{px +uf +Py^ -ml , 
which can never be fulfilled. The second solution, i.e. u=px, results in 



Pv 



(70) 



(71) 



(72) 



(73) 



Therefore, no poles shall occur for \v\ < m2, while for \v\ > m2, the integral (68) becomes ill-defined. This is important 
as it shows us that also in the well-studied case of two real masses, the original integral (68) is also only well-defined in 



a certain region of the complex plane, here displayed in FIG. 5. Consequently, one also needs to perform an analytic 
continuation outside this region, just as in the case of pure complex masses. 




FIG. 5: Parabolic region in the complex plane where the momentum integral 02{k^) exists with ml = 1/2. 



HIGHER LOOP CASE 



Having discussed in detail the analyticity properties of one loop correlation functions made of i-particles, we face 
now the extension to higher orders. As we shall see, the use of i-particles turns out to be extremely useful here, 
enabling us to construct a rather interesting iterative procedure valid for an n-loop integral. Note that in the toy 
model only one type of diagrams appear, sometimes known as water melon diagrams, because there are no interactions. 
The case of diagrams including interactions as appearing for the GZ theory is left for future studies. Let us begin 
with the two loop order. 
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A. A two loop example 

In order to obtain an example of a two loop correlation function, we consider the local operator 

02(x) = X{x)7j{x)U{x) . (74) 
For the correlation function {02{k)02{—k)) we find 

(02(fc)02(-fc)) = J2{k^) = / ttt^^t;^^ 7= ^ , (75) 

^ ^ ^ ^ ^ ' J (27r)2 (2^)2 g2 p2 ^ ,^202 (fc _ g _ p)2 _ ,^202 ' ^ > 

where k?' is the Euclidean external momentum. In order to evaluate it®, we rewrite it as 

'^^^^^^ " J J2^f (/ (27r)2 p2 + ,^02 (fc _ ^ _ p)2 _ ,^202 ) ' 



so that we can use the spectral representation ( 45 ) already found at one loop level, namely 

dr p{t)—-. , (77) 



(27r)2 p2 + ,^02 (fc _ g _ p)2 _ j^02 '^^ V + (fc - g)2 ' 

with To = 2726*2. Therefore 

1^^,^^). (78) 

The q-integral becomes now straightforward. It corresponds to the one loop integral in which one particle is massless, 
and the other one has a real and positive mass r. Such integrals also have a spectral representation with positive 
spectral density, given by 

y (2^?r+(fc-,)2=X (^9) 
where pi{s) is the associated spectral function. Its explicit value can be found in Appendix [A| eq. (B6). Thus, 

/oo poo 1 

dTp{T)j dsp,{s)-^. (80) 

This is not yet in the form of a spectral representation. By switching the order of integration, we can however reexpress 
the double integral as 



oo 



Mk')^ I ds^— r drp,{s)p{r)^ r M±ds. (81) 



TO s + k'^J^g y.ro s + k^ 



This means that the spectral density is given by 

P2(s)= / dTp(r)pi(s). (82) 



whereby we kept in mind that pi{s) will depend on t as well. This allows to make the following observations: 

• The branch cut will be located on the negative real axis from — oo to — tq. 

• For s £ [tq, oo], we have p{t) > due to the positivity of p(r) over the interval [tq, oo]. We shall thus have that 
P2{s) > itself if pi{s) > for s E [to,oo]. The latter turns out to be true, as pi{s) is defined for all r with 
T > To, where it is positive. 



Prom now on, we shall not bother about potential subtractions to make it well-defined. These can always be obtained by taking a 
suitable number of derivatives w.r.t. the external momentum k^. 
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In conclusion, the two loop correlation function (75 1 can also be reexpressed in a spectral form, with associated 



positive spectral density. As we have seen, there is even no need to explicitly evaluate the integral or spectral density 
to establish these conclusions. To completely convince ourselves, we can also evaluate (82) explicitly, by making use 
of (46) and (B6), which leads to 



P2{s) = 



Vt2 - 8614 




2V2e'^ ^arctan 



86*4 




(83) 



and then one can check that p2{s) is indeed positive for s > 2^/29^ 



B. A three loop example 

To construct another example of a higher loop correlation function, one might consider the composite operator 

Os{x) = \iX{xMx)f , (84) 

which gives rise to a three loop expression 

Jffc2)=/^^^^ ^ ^ ^ (85) 

' J i2nY {2ny {2T:)d q2 i^02 p _ ^^02 p2 ^^02 ^ q _ I _ p)2 _ ^^29^ ' ^ ' 



We proceed in a similar way as for the two loop case. Starting from 



and by making the change of variables 

e^£-q, (87) 



we get 



Thus, 



from which 



J3{k')= / drpir 



d'^e 1 f d'^q 



{2t:Y T + {k-lY J {2TT)d q2 + ^72612 (£ - qf - i,/29^ 



oo 



Mk')= I drpir) I ,f t , ,1 r dsp{s)^ , (89) 



To 



POO POO POO -f 

Me)= drpir) dsp{s) drpi(r)— — (90) 

follows, keeping in mind that the branch cut for two particles with respective positive real masses p^ and P2 starts 
from —{pi + P2Y; see also Section IV B 



We see that we again find a kind of convolution of already known spectral representations. As before, we 
still need to show that J^ik^) can be brought into a spectral representation with positive density P3(r) for r > r* > 0, 
thereby establishing a physical branch cut along [—00, — t*], with yet unknown. It turns out that this is little more 
complicated than in the two loop case. Let us look at the underlined part of (90) first, which is a generalization of 
the previous case J2{k^)- Hence, we wish to examine 

/•oo I'OO 

Mk^)^ dsp{s) drp,{r)—-, (91) 
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where a > ^Jtq- We propose the substitution = \fs + a, leading to 



where we have introduced a novel spectral density 



dim \ drp^{r) 



r + fc2 



^ri- 



(92) 



(93) 



It is clear that for ^ > + a)^, we have p(^) > 0, invoking the properties of p(s). We then observe that we have 
managed to reduce J2(fc^) to the J2{k^) case already studied, see eq. (801, meaning that (91 1 has a spectral density 
representation of the form 



1 



As such, we can write 



Mk^) 



with P2{s) > . 



This is again of the type J2{k'^), from which we conclude that we can write 



Uk^) 



with p3(s) > 0, for s > t* = 8\/26'2. 



dsp-i{s)- 



1 



fc2 



dspsi-s)- 



4to 



(94) 



(95) 



(96) 



We may summarize by noticing that the higher loop correlation functions can be written as a kind of convo- 
lution of the lower loop spectral representations. If the involved lowest, i.e. two loop, spectral densities are physical, 
so will the higher loop spectral densities be, by making use of the iterative argument given in this section. 



VI. A FEW REMARKS ON THE PHYSICAL MEANING OF THE i-FIELDS (A, -q) 

As we have seen in the previous sections, the introduction of the «-fields (A, rj) enables us to construct in a relatively 
simple way a set of composite operators whose correlation functions exhibit real cuts only. It is worth thus spending 
a few words on the physical meaning of these fields. 



Looking at the very starting point, eq. (18 1, we see that the action of the field V has been modified by the 
addition of a nonlocal term, accounting for a deep modification of the long range behavior of the model. It is easy 
to figure out that the field ip itself acquires thus a long range nonlocal component. The construction of a suitable 
set of composite operators in terms of the field ip becomes a highly nontrivial task, as nonlocal terms have to be 
incorporated in order to achieve meaningful correlation functions. When the model is cast in local form, the long 
range behavior of the theory is accounted for by the auxiliary localizing fields ip, (p. The fields ■(/;, Lp, Cp have to be 
considered thus on an equal footing. It is instructive now to express the composite operator Oi{x) = X{x)'q{x) in 
terms of the original field tj}. From the equations of motion, it follows 

9^ 

Analogously, for the i-fields, one gets 
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1 / (P' 

Therefore, for the operator 0\(x) one obtains 



(98) 




0^{x) = X{xMx) = + 2 (^(fg^V^ ) j . (99) 

from which one clearly sees that nonlocal terms are needed to achieve a sensible operator. The i-fields A, 77 provide 
thus the correct field variables which enable us to construct in a simple way, and within a local quantum field theory 
framework, the relevant composite operators of the theory by taking into account the nonlocal long range effects. 



VII. INTRODUCING i-PARTICLES FOR THE GRIBOV-ZWANZIGER ACTION 

Let us discuss here how i-particles can arise in the Gribov-Zwanziger action, given by expression (|8|. In what 
follows, we shall limit ourselves to evaluate correlation functions of suitable composite operators at one loop order 
only. To introduce the i-particles in the Gribov-Zwanziger action, it suffices thus to consider the quadratic part of 
expression (|8| containing only the fields (p"'', (^^''), namely 



where use has been made of the transversality of the gauge field, 9^^^° — 0. We now proceed by decomposing the 
fields i'Pa^j'Pn^) in symmetric and anti-symmetric components in color space, i.e. 



vi^'^^U^f+v';) , (101) 



and 



Thus 

A first step towards diagonalization of this expression is achieved by setting 
so that 

(105) 
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From expression (1051 one sees that the gauge field mixes with the adjoint projection of Vjl"''', obtained by 
employing the following decomposition 

TA[a?j] rabp rprnn-tr[rnn] , I ^r[ah] rahp rpra7i^r[mn] \ f^^^PT/P i cff*^] 



where 



-trp _|_ rpmn-tr[7nn] 



(106) 
(107) 



stands for the adjoint projection of V^l"''' in color space, and 

g[ab] _ ylab] ^ jabp jpmnylmn] (108) 

denote the remaining independent components of V^"''' which are orthogonal to the tensor f"-^'^. In fact, making use 
of 

jahc jdbc j\^^ad 



it is easily checked that 



fabcg[ab] ^ q 



Therefore, expression (105) becomes 



Finally, setting 



one obtains the diagonal action 



The fields A° , 77° describe the i-particles of the Gribov-Zwanziger action. 



(109) 
(110) 



(111) 



(112) 



(113) 



A. Evaluating correlation functions in the Gribov-Zwanziger theory using i-particles 

Let us give here two examples of one loop correlation functions of composite operators constructed from the i-fields. 
To this order, one can introduce the i-field strengths defined by 



(114) 



The propagators are^ 



{K{k)\l{-k)) 



^ab 



There are more propagators than the ones shown, but these are not relevant for the calculation presented. 
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~fc2~ 



with 7** = 2 7^ as introduced in Section |lj 

As simplest examples we investigate the following composite operators at leading order: 



(115) 



O 



The integrals to be calculated for the correlation functions are 

{0-mO-H~k)) = iN^-l) ' ^V\v-kf+^id-2)i,^-,kf 



(0^^>(k)0'^>(-k)) (N' 1) I 32(fcV_(fcp)2) 



(116) 

(117) 
(118) 



To avoid to dwell on technical details, let us here only outline the results for d = 2. We have collected the technical 
details, including those for d = 4, in Appendix [C] We are interested in 



i(i) 



with 



F{k') = 



d'^p p^{p-kf + {d-2){p^ ~pkf 
{2ttY {p^ + ij^){{p - ky - ij^) 



(119) 



(120) 



After the introduction of a Feynman parameter and some manipulations, it turns out that we may write, see (C17l, 

"9F(P) 



F{k') - k' 



1 

An 



+•00 



dT 



dk^ 
1 



1 



2^/7 



1 fc2 



We thus find 



-m 

T 47r 



Pir) 



+00 



dT 



1 



1 



2V7 



1t2 



1 

in 



+00 



dT 



1 



1 



2^1 



1 T 



We conclude that, upon restoring units, we formally have 

(oS(fc)o5(-fc)> = 



+00 2(^2 _ ^^^4 
2 ny/T^ -ij^T + k^ 



(121) 



(122) 



(123) 



2f 



We clearly notice that the spectral density pir) is positive for r > 27 . The result as written in (1231 is indeed 
only formally correct, since the l.h.s. of (1231 is divergent, directly seen upon inspection of its definition (1201. 
Nevertheless, the spectral representation ( |123 1 appearing in the r.h.s. defines a finite function. The apparent 
contradiction is easily resolved by realizing that one should in fact refer to ( 121 1, which gives the correctly subtracted 
result. 



Let us now turn to the analysis of 



(0S(fc)0rj(-fc)>=32(iV^-l)G(fc^), 



1(2) 



with 



G(fc2) = 



d'^p k'^p^ - [kpf 

{2nY {p^ + ij^){{p - ky - ij^) 



(124) 



(125) 
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We again refer to Appendix [C| After the smoke clears, we find in (C28l 



G{k') - G 

\ ^+00 



The spectral density can be read off, 



-G(0) 

1 '•+°° 



.+00 ^ ^ .+00 

/ dry/r^ - 1— / drv/T2 - 1- 

Ji Stt 7i t 



P{t) = ^Vt2 - 1 , 

whereby p(t) > for r > 1. We reintroduce the units, and we conclude that 



(oS(fc)og(-fc)) 



4(iV2 - 1) 



+00 



dr 



272 



- 

T + fc2 



(126) 



(127) 



(128) 



which is again a formal result due to the divergent nature of both l.h.s. and r.h.s. . 

Using a similar analysis, one can also derive the spectral densities if d = 4. The spectral densities in all 



cases turn out to be positive and thus at least at leading order the operators O^j* and O]^' appear to be physical 



1(2) 



Ar, 



'Ar, 



We conclude by remarking that the physical meaning of the i-fields A",?]" is akin to that of the i-fields (A,?]) 
of the toy model. In fact, from equations (111 I, (112) one easily gets 



A'^ = — 



'2N 



(-92) 



(129) 



and a similar expression for ry^j,. As expected, the i-field strength AJJ^ contains an explicit dependence from the 
nonperturbative Gribov parameter 7. This dependence signals that, in order to construct a sensible operator displaying 
good analyticity properties in the Gribov-Zwanziger theory, the starting gauge invariant operator obtained in the 
Faddeev-Popov theory, i.e. without implementing the restriction to the Gribov region 17, has to be supplemented 
by the addition of terms which are 7-dependent. To some extent, one might figure out that a would-be physical 
operator in the presence of the Gribov horizon has to be constructed by deforming the starting gauge invariant 
operator by appropriate terms which exhibit an explicit dependence from the Gribov parameter t, and which allow 



for a cancelation mechanism of the unphysical cuts appearing in expressions like those of eqs.(14), (15). Although the 



I- variables (AJ^^, 77^,^), introduced here in the quadratic approximation, enable us to obtain in a relatively simple way 
examples of correlation functions possessing a Kallen-Lehmann spectral representation at one loop order, it is certainly 
tempting to look at those variables in the full Gribov-Zwanziger action, trying to unravel a systematic mechanism to 
take advantage of the presence of the Gribov horizon to cancel the unwanted unphysical cuts. 



VIII. CONCLUSIONS 



In this work we have pursued the investigation of the analyticity properties of correlation functions evaluated 
with a confining propagator of the Gribov type. In particular, as illustrated in the toy model, we have been able to 
characterize examples of composite operators whose correlation functions display cuts only on the negative real axis, 
while possessing a positive spectral function, a result which has been extended to higher loop correlation functions. 
The introduction of i-particles, which seem rather natural objects when dealing with a Gribov type propagator, has 
proven to be very useful in the construction of such composite operators. 

In the case of the Gribov-Zwanziger theory, so far, we have been able to provide examples of composite oper- 
ators, made of z-particles, whose correlation functions at one loop order exhibit real cuts only. 

The introduction of i-fields in the full Gribov-Zwanziger action within a local and renormalizable framework 
is certainly a point worth to be investigated. This could open the possibility to obtain examples of correlation 
functions displaying good analyticity properties at higher orders, a result which can be regarded as a highly nontrivial 
achievement. The systematic characterization of how a given gauge invariant operator needs to be deformed by the 
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Gribov horizon in order to possess a Kallen-Lehmann representation is a big challenge, requiring many ingredients as, 
for example, the mastering of the renormalization procedure of gauge invariant operators within the Gribov-Zwanziger 
action, a necessary step for a consistent higher order calculation. Let us remind here that, due to the soft breaking 
of the BRST symmetry, the issue of the renormalization of gauge invariant composite operators requires a careful 
analysis, as recently done in [26 , where a renormalization group invariant, and thus a fortiori finite operator in the 
Gribov-Zwanziger theory, which contains the scalar glueball operator F'^{x) — F°^(a;)_F'^y (a;), was already identified. 
In particular, we point out that the soft breaking of the BRST symmetry implies that the operator mixes not only 
with BRST exact composite operators, but also with BRST non-invariant local quantities, which are determined by 
the softly broken Slavnov- Taylor identities |26| . Evidently, the situation generalizes to higher dimensional local com- 
posite operators containing three or more field strengths F^^^ which are contracted in such a way to give rise to a color 
singlet as, e.g. , /"''"^Fj^^^Fl^pFp^. Besides the mixing with other gauge invariant operators, a nontrivial mixing matrix 
between BRST exact terms as well as BRST non-invariant terms is to be expected here too. This mixing should be 
taken into proper account, as (1) it is indispensable to obtain finite results, (2) it can play a major role in identifying 
the analyticity properties of the complete operator, which differs from the classically gauge invariant one. From this 
perspective, the lowest order results presented for the Gribov-Zwanziger correlation functions (C41 1, (C42l, (C43l and 
(C44l are only a first step. It will be interesting to find out whether one could also construct a renormalization group 
invariant operator in terms of the i-particles, and what the connection would be with the one already discussed in |26] . 



As this paper only contains the first effort in constructing meaningful correlation functions in the presence of 



a Gribov type propagator (13 1, we did not yet study the more general gluon propagator 



{K{k)AU-k)) = 5- (S,. - ^) , = , (130) 

which was discussed in |12l UM giving rise to the so called refined Gribov-Zwanziger action. It arises when 
additional nonperturbative quantum effects are considered. One might wonder if it would still be possible to obtain 
correlation functions with only physical cuts along the negative real axis when referring to the gluon propagator 



(130 1. One can easily write down the toy model analogy of the refined Gribov-Zwanziger action. It turns out that 
the previous question can be answered positively, again by introducing the i-particles. We refer to [30:, as the explicit 
evaluation of the spectral densities calls for a set of mathematical tools rather different from the ones used in the 
current paper. 



Let us end this work by mentioning that several aspects related to the Gribov-Zwanziger theory as well as to 
the use of the Gribov type propagator and of the corresponding i-particles interpretation are still open. This is the 
case, for example, for the construction of the Minkowskian version of the theory, which remains to be achieved. We 
emphasize that, as it stands, the Gribov-Zwanziger action has to be understood as an Euclidean field theory whose 
origin can be traced back to the lattice formulation of gauge theory. 

It is worth recalling here that gauge theory is well-defined on a finite lattice. It is an Euclidean theory for 
which the correlation functions of gauge non-invariant fields such as the gluon and quark fields vanish [31 . The 
Osterwalder-Schrader refiection positivity |17| [TH] holds for gauge invariant operators, so that for these quantities 
there exists a positive metric quantum mechanical Hilbert space [1311 132|. Moreover, there exists a positive transfer 
matrix T, with < T < 1, which in the continuum limit would imply a positive Haniiltonian |31fl33| . In |34| . it has 
been argued that the critical limit of lattice gauge theory, in the minimal Landau gauge, is precisely the theory we 
are discussing. If so, the Gribov-Zwanziger theory would inherit from lattice gauge theory a positive metric Hilbert 
space for composite gauge invariant quantities and a positive Hamiltonian, although in four dimensions this is not 
trivial because of renormalization. 



We also point out that the complex singularities exhibited by the Gribov propagator and by the i-particles 
might jeopardize the usual implementation of the Wick rotation in the correlation functions, so that the construction 
of the Minkowskian version of the theory is still an open point. In particular, as is apparent from the presence 
of complex poles, the two-point correlation function of the elementary Euclidean gluon field, as described by the 
i-particles, cannot be interpreted as a propagator of a physical particle in Minkowskian space, but may be appropriate 
for a confined gluon. On physical grounds, we expect that the rotation to Minkowskian space should be possible 
only for a restricted class of composite operators, which should be related to the physical spectrum of QCD. To 
some extent, this is precisely what emerges form the present analysis, namely: the Kallen-Lehmann representation 
seems to exist only for a restricted class of operators. The lack of a Minkowskian description of the i-particles 
and the presence of complex poles also indicate that the usual assumptions (or axioms) of local quantum field 
theory do not apply, in particular the spectral condition about the location of branch cuts for a two-point function [35] . 
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The Euclidean nature of the i-particles allows us to argue that they describe a phase for which the gluons 
are unphysical, i.e. they do not appear as asymptotic states in the 5'-matrix elements. As such, the issue of the 
cancelation mechanism of the singularities in the S-matrix giuon amplitudes |36H40| . which would arise due to the 
use of a propagator as (13) or (130), is not addressed here. Evidently, we are still far from a concrete description, 
within a purely quantum field theory framework, of the mechanism which would account for the conversion of the 
gluon cloud into physical jets of hadronic matter^''. 



In conclusion, although the results which we have obtained so far look promising, a satisfactory characteriza- 
tion of the analyticity properties of the correlation functions of the physical composite operators which would 
correspond to the spectrum of QCD remains to be achieved. 
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Appendix A: Kallen-Lehmann representation of the toy model two-point function {Oi{k)Oi{—k)) for general d 



1. Using the Schwinger parametrization 



In this Appendix, we shall verify that the two-point function {Oi{k)Oi{—k)) , defined in eq. (32), has a spectral 
representation with associated positive spectral density for any value of d = 2, 3, 4. We have already discussed the 
case d = 4 in section [III B| but we shall follow a follow a different route here, which is applicable for general d, thereby 
serving as a check on the d = 4 result we derived in expression ( 46 1 . Returning to expression ( 32 1 , we observe that 
by using the Schwinger parametrization, given in eqs. ( |C29 ) and (C30l, it can be cast into the form 



(Oi(fc)Oi(-fc)> 



d'^p p^{k-pf + 2e^ 



(27r)^ ((fc-p)4 + 26'4)(p4 + 26l4) 

da / dp cos(V29^{a- (3) 
Jo ^ 



oo /"OO 

da / df3 cos(V29^{a- (3) 
Jo ^ 



1,2 

e 



da / dp Re' 
Jo 



(47r(a + /3))- 
(47r(a-f J ' 



(Al) 



The cut structure is encoded in the exponential of the last line. In order to make this structure more apparent, a 
few algebraic manipulations of this expression are needed. Following as a first step we insert in (Al) the unity 
decomposition 



dXS{a + P-X) 



(A2) 



See also the related discussion about the possible interpretation of the gluon dispersion relation implied by the use of a Gribov type 
propagator presented in |41| . 
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After a rescaling a — ^ Aa, (3 — ^ A/3; the integral over A can be evaluated, yielding 
(C'i(fc)Oi(-fc)) = / da / d/36{a + (3-l) I dX X^-^^Rc ' 

da 



r(2-f) 



(47r(a + /3))5 J 
dl3S{a + l3-l)ReUa(3k'^-iV20^{a-(3)y ^| . (A3) 



(47r)2 

The second step is to parameterize the branch cut. For that purpose, exactly as before, we introduce the variable u 

_ a - /3 2a - 1 



2a/3 2a(l-a) 



(A4) 



where the last equality follows from the constraint imposed by the delta function: /3 = 1 — a. In terms of u, expression 
(A3) takes the form 

(Oi(fc)Oi(-fc)) = ^ /^ Re / du— =(-——=] (k'-t2V2e'u)' . (A5) 

Again we note that due to the first square-root in the integrand, this expression has a branch cut in the complex 
M-plane which starts at m = i and extends till u = ioo. This cut may be rotated to the real axis by exploiting the 
properties of the contour integral to express the result as an integral along the cut exactly as we did before. Defining 
M = it/ we get 

(0,(fc)0,(,fc)) ^ ^^^^:J^ Rel Tdy-^^^ ( ^^_ \ ^ (fc^ + 2^/20^y) " ^ . (A6) 



For the case d = 2, we can infer the spectral density directly from (A6I. We find 



hence 



1 dr 1 

{0,{k)Oi{-k)} = - / , (A7) 



with indeed Pd=2{T) > for r > 2^/26^. Explicitly, the two-point function yields 

2ttV80* - k* 



(Oi(fc)Oi(-fc)) = ^ : , (A9) 



displaying a branch cut from —00 to —2\/2(P' . 

For a general d, it is useful to employ the following identity 



^-2 1 (T-2y2f?2y)^ ' 



which already displays the familiar structure found in the Kallen-Lehmann representation. Inserting this identity in 
( A6 ) , it is not difficult to see that the resulting overall range of integration can be reordered as 



We make now a further change of variables defining 



dy / dT(...)= / dT / '^'"dy(---) . (All) 



1 TT 

- = C0S7; < 7 < - . (A12) 

y 2 
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Putting all together, it is just a matter of algebraic manipulation to obtain 



(Oi(fc)Oi(_fc)) = r dr^^ , (A13) 



where the spectral function p(r) is given by 



p{t) = — — ^ / ^7 cos ( (^ — 1)7) (cos 7 — cos (/"(r)) ^ , (A14) 

with 0(t) defined as 



cos0(t) = ?^. (A15) 

T 



For d = 4, the spectral function can be immediately evaluated 



Pd=4(T) = -^sin0(r) - -i-2Wl- ^ , (A16) 



which is consistent with ( 46 1 . 



When d — 3, we need to evaluate 



^ ^ 1 1/"^^ cos(7/2) , 

4V27r2 Jo V cos 7 - cos 



We must compute an integral of the class 



cos2"+i(7/2)d7 _ 2 Z^^' cos^"+i(y)dy ^ 2 /'"^^^'^ (1 - f3^rdl3 ^^^^^ 



Vcos7-cos(/) Jo - 2 sin^ 7' - cos ^ y^l- cos 0-2/3 



with n e N, where we used the substitution sin 7' = /3. This integral ( A18 1 can be expressed in terms of the Gaussian 
hypergeometric function, 

^ - 2^^1(1/2; -n;l,sin2(0/2)). (A19) 



Vl-cos<^-2^2 ^ 
Since, in particular 2^1 (1/2; 0; l,x) = 1, we obtain 



P<i=3(T) = ^^, (A20) 



also positive for r > 2^/26^. As such. 



F(.^)./ _ ^r^ = i4arctanJ^, (A21) 



72^262 T + fc2 47r fc2 y fc2 

for d = 3, valid in the complex plane, exhibiting a branch cut for k"^ G [^00, —2y/29 



2. Using hypergeometric functions 

A very convenient method to evaluate Feynman diagrams is the use of hypergeometric functions, since these 
are meromorphic functions in the complex plane. The matter of analytic continuation is especially clear, when 
one uses these functions as will become evident in the following calculation. Another advantage is the possibility 
of keeping more parameters general like the number of dimensions or the exponent of the propagator, which can 
even be non-integer as needed sometimes. Hypergeometric functions have proven useful especially for higher n- 
point functions or diagrams with many mass scales, because they allow to write down the result as a closed expression. 
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One possibility to arrive at hypergeometric functions is the use of the negative dimension integration method 
(NDIM) |42fH5] . but one can for example also work with Mellin-Barnes representations |3B] or dispersion relations 
|47| . The result for the one loop two-point diagram with masses^^ := rnf = i\f2(P' and = ~i^/29^ in d 
dimensions with exponents ii and Z2 for the propagators is |48j 



\ (-"/2 -l2~ l2,d/2 + 12) 

X Fi{-n -12- d/2, -^2; 1 - Z2 - d/2, d/2; — ) 

imlf/'+'- {-12, -d/2)Fi{-i2,d/2; l + d/2 + i2,d/2; 4> (A22) 

7712 m| y 



where the Pochhammer symbol is defined as 



T(a + n) , . ^ 

ia,n)^ ^ (A23) 
r(a) 



F4 is the fourth Appell function, which is a meroniorphic function of two variables. At this point we will choose a 
specific series representation for it: 

^, , , , ^ {a,m + n){b,m + n)x"^y"' /^r,^^ 



This restricts the validity of the following manipulations to a certain region of convergence, namely y'lx] + y'\y\ < 1. 
The fact that our masses are complex fits quite naturally in this framework, as the variables can be complex 
numbers. Later on we will switch back to hypergeometric functions, where we do not have to worry about regions of 
convergence, since we can always choose a series representation appropriate for the values of the variables x and y. 

We continue by rewriting the Appell function into a Gaussian hypergeometric series, 

Fiia, h; c, d; a;, y) = V ^"'"^^^.'"^ ^ 2^^1(0 + n,h + n;d; y), (A25) 



and evaluating the one-dimensional series with 



We arrive at 



,F^{a, b;a+l-b;-l) = 2-°^i/^ , , | , 

^ ' ' ' ' r(l/2 + a/2)r(l + a/2 - 6) 



/2(-l,-l;mi,e'"/2mi;d) 

_U ^-d/2 l/2( / 2Nd/2-22-2+d/2 {d/2, -l)T{2 - d/2.) 

^ \^ ^' (-d/2 + 2,d/2-l)r(3/2-d/4)r(l-d/4) 



(A26) 



E 



(-d/2 + 2,n)(l,n) 1 ( -k^\ I 



{d/2,n) {i/2-d/A,n/2){l-d/A,-n/2)\2m? J n\ 



+ {.r?rH.r?)''^~'2'''^ {l,-d/2)T{d/2) 

+ ^""1^ ^""2^ ^ r(i/2 + d/4)r(d/4) 



The choice of the square root of two for the masses is in correspondence to the main part of the paper. 
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^ {l,n){d/2,n) 1 fz^Yl.] 

^ (d/2,n) (l/2 + d/4,n/2)(d/4,-n/2) I 2m2 i n\ ' ^ ' 



oo 

X 

n=0 



id/2,n) (l/2 + d/4,n/2)(d/4,-n/2) \2m^ J 
To deal with the half-integer summation indices, we split the sum into even and odd indices: 
72(-l,-l;mi,e''^/2^i;rf) 

=(47r)-'^/27ri/2 ( (rn2)'^/2-22-2+<^/2- '^''^^ "^^^^^ " ''/'"^ 



(-d/2 + 2, d/2 - l)r(3/2 - d/4)r(l - d/4) 

/y^ (-d/2 + 2,2n)(l,2n) 1 f -k^V" 1 

^ (d/2,2n) (3/2-d/4,n)(l-(i/4,-n) V2m2j (2n)! 

(-d/2 + 2,2n + l)(l,2n+l) 1 /-fc2\ 2"+^ 1 

(d/2,2n + 1) (3/2 - d/4, n + 1/2)(1 - d/4, -n - 1/2) 1^2m2 j (2n+ 1)! 

I . 2^-1. 2^./2-i.-./2l (l,-d/2)r(d/2) 
+ (mij (m^j 2 4r(l/2 + d/4)r(d/4) 

/y^ (l,2n)(rf/2,2n) 1 /-PN"" 1 

(d/2,2n) (l/2 + d/4,n)(rf/4,-n) V2m2j (2n)! 



X 

n=0 



(l,2n + l)(rf/2,2n + l) 1 /-py^+i 1 

"""^^ (d/2,2n+ 1) (l/2 + rf/4,n + l/2)(d/4,-n- 1/2) V2m2j (2n + l)!*^ 

-r47rr''/27ri/2 ( /™2Nd/2-22-2+d/2 {d/2, -1)T{2 - d/2) 

' y ^' (-(i/2 + 2,d/2-l)r(3/2-d/4)r(l-d/4) 

X (2^1(1, l-rf/4; 1/2 + d/4;-^) 

^-^/^"^'^^ .^.(l,3/2-./4;l../4;-^)) 



2m2 (rf/2,l)(3/2-d/4,l/2)(l-d/4,-l/2) ^ ' ' / ' ' / ' 4^4^ 
+ (-D-Hm^)^/-^2-/2^(i^^^^^ ,^,(1, 1 - ./4; 1/2 + ./4; -^) 

2Fi(l,3/2-rf/4;l + d/4;-^))j, (A28) 



2rn2 (1/2 + rf/4, l/2)(d/4, -1/2) 
where the following identities for Pochhammer symbols were used: 

{a,m + n) —{a,m){a + 111,11,), (A29) 
{a,2n) =22"(a/2,n)(l/2 + a/2,n), (A30) 
r{a + m) =r{a){a,m), (A31) 
(-1)-" 

{a,—n)=- ^, ninteger. (A32) 

(1 — a, n) 

For specific values of d the expression will simplify. 

Here wc arc dealing again with a hypergeomctric function. Wc know that 2Fi{a,b;c;x) is analytic in the 
complex plane except at +1. from where a branch cut to +00 starts^2 From the argument of 2.P1, it looks as if the 
result has two branch cuts: From ±.2^/26"^ to ±00. We will see below in the expressions for specific dimensions that 
only one of those remains, namely the physical one on the negative real axis. 



Strictly speaking this is only true for nonexceptional parameters o, h and c; these are parameters for which the series terminates or the 
series is undefined. 
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Now we want to obtain specific expressions in two, three and four dimensions. In two dimensions we get 
/2(-l,-l;mi,e*''/2mi;2-2e) 



(47r)-i+^7ri/2 I (mi) 



X (2^1(1, 1/2 + £/2;l-e/2;-^) 



(l-e,-l)r(l + £) 



(1 



2m2 (l-£,l)(l+e/2,l/2)(e/2,-l/2) 
^ r(l-e/2)r(l/2-e/2) 



(1 + £, 1 - £)r(i + £/2)r(i/2 + £/2) 

2i^i(l,l + e/2;3/2-£/2;-A_ 



))- 



fc2 



2m2 (1 - £/2, l/2)(l/2 - £/2, -1/2) 



(2^1(1, 1/2 + £/2;l-£/2; 



2^^1(1, l + £/2; 3/2- £/2;- 



(A33) 



Although it is not evident from above due to divergent coefficients, the integral is finite in two dimensions. Performing 
the hmit e — J' yields 



/2(-l,-l;TOi,e*"/2mi;2) 



7T/2 - arcsin (^) ^ arccos (^) 



27rV86i4 - 



(A34) 



The arccos has branch cuts from ±1 to ±00 and the inverse square root from to —00. They combine in such a way 
that only a cut on the negative real axis starting at —2\^9^ remains. 



For three dimensions the limit £ — !■ yields 



/2(-l,-l;TOi,e*"/'mi;3) 



6-\/26'2 2-^1(2, 1, |, ^) — 2-Fi(|, 1, |, ^] 



487r6'3 

With the integral representation for the Gaussian hyper geometric series, 

r(c) 



2Fi{a,b, c,x) 



r(6)r(c-6) 

this can be rewritten into a spectral representation 



/2(-l,-l;mi,e*^/2mi;3) = 

87r 



f dtt''-^{l~tY-^-^{l~tx)-'', Rec>Re6>0, 

JO 



(A35) 



(A36) 



dr- 



.-1/2 arctan( 



fc3 
2%/2e2 



2^292 



r + A:2 



47rvfc2 



(A37) 



As the inverse tangent has branch cuts from ±z to ±joo, we again observe a branch cut on the negative real axis. 
The role of the square root in the denominator is to compensate the cut stemming from the square root inside the 
arctan, which would start at fc^ = 0. We notice the expression (A37l coincides with (A21l. 



Finally we evaluate eq. (A28l in four dimensions. Here more care is necessary since the integral is divergent: 



/2(-l,-l;mi,e'"/2^i;4-2e) 



r(i-e)r(e)r(e) 



2^-.oel r(e-i)r(2-e) 



4r(3/2-e/2)r(l-e/2) 



X 2Fi(e/2,l;3/2-./2;- — ) 



r(2-e)r(l/2 + e/2)r(e/2) 

Am- 

2.-enA( r(e-l)r(2-e) r(6 + l)r(2-e)r(l/2 + 6/2)r(6/2) 
2to2^ 4Vr(3/2-e/2)r(l-e/2)r(e)r(3-e)r(l + e/2)r(-l/2 + e/2) 

r(6-l)r(2-e) r(l-e/2)r(3/2-6/2) x 
r(3/2 - e/2)r(l - e/2) r(2 - e/2)r(l/2 - e/2)J 
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X 2Fi(l,l/2 + e/2;2-e/2;- — ) j. (A38) 

The coefficient of the second hypergeometric function and the function itself are finite for e — >■ and only the first part 
yields divergences. To extract one factor of e we can use the integral representation of the Gaussian hypergeometric 
function, eq. (A36l. For the given values of the parameters (a = 1, 6 = e/2, c = 3/2 — e/2) there is a pole at i = in 
the integrand. Rewriting the integral as 



(A39) 



with 



J^ = ^J}f^ M [' dtt'-\l-tr-'-^f{0), (A40) 



mnc-b) Jo 

it can be exposed explicitly. For the last term we calculate 



r(6)r(c-5) Jo 

^2 = ^J^ f dtt'-\i-tr-'-\f{t)-fm, (A41) 



f{t)~f{Q)=t——^tg{t). (A42) 
i — tx 



This raises the exponent of t in J2 by one. The integral in Ji is identified as the /3- function, 



Jo r(2; + y) 



(A43) 



Thus it yields Ji — 1. Note that the conditions for this identity coincide with those of the integral representation of 
2F1. The second integral is written back into a Gaussian hypergeometric function: 

mnc-b) Jo ^"(1-^) T37^-^r(6)r(c-6) — r(^TT^ — 2F,{i,b + i,c + i,x) 

= ^^x,F^il,b+l;c-b;x) ^ - -^^Z^ ,^^,(1, e/2 + 1; 5/2 - e/2; -^). (A44) 



The coefficient of this function in eq. (A38l has a term of order 1/e. Thus we need the total expression in eq. (A44l 
up to order e. The goal in rewriting the original Gaussian hypergeometric function like this was to extract one factor 
of e, because now we need only the lowest order in e of the Gaussian function. Combining all intermediary results we 
obtain 

, , ,,/2 , 11 2-7 + ln(4^)-ln\/2e2 ^^,^02 + ^gQ4 _ ^4 arccos(fc2/(2V202)) 
J2(-l,-l;mi,e / m^; 4 - 2.) = ^ - + + . 

(A45) 



Again the cut the on the positive real axis is canceled by the one from the square root. The subtracted result is 

167r2fc2 



w 1 1 ^n/2 , o^ fc^ - nV29^ + _ ^4 arccos(fcV(2V202)) 
J2(-l,-l;mi,e*^/^mi;4-2e) = — — ^ '-^ ^, (A46) 



in accordance with (47 1. 

The explicit example of the one loop case demonstrates that hypergeometric functions are a useful method, 
which can develop its full power especially for higher loop cases. The derivation for an A^-loop diagram in terms 
of generalized Lauricella functions can be found in [46^ . For diagrams with only two massive propagators only two 
hypergeometric functions remain that can be treated along the lines above. The end result is again a cut only on the 
negative real axis. 



Appendix B: Spectral density for one real mass and one vanishing mciss 



We start from 



{k — + rrfi ' 

For fc^ > 0, we can employ the Feynman trick and differentiating w.r.t. k"^ yields for d = 4 



dF{e) 



dx 



x{l — x) 



x{l — x)k'^ + xm? 



dx- 



167r2 Jo k^ + f: 

2 

We perform a transformation of variables, by setting s = f^, 



dF{k^) 

After doing a partial integration, we obtain 

dF{e) 1 



167r2 



+ 00 



ds- 



1 1 



fc2 
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'""^ Jm'^ ds \ s J k^ + s ' 



dk'^ 



167r2 

-1 
167r2 



1 m 
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fc2 + S S 
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A;2 + 
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ds 



I ds— — ( 
J ^2 s ds \k^ + s 



1 



s {k'^ + sf 



-1 d 



SO that 



F{k'^) - F{0) 



1 

167r2 
1 

167r2 



-In(fc2+TO2)+Inm2- / 

J rr 

r+oo ^ p+oo ^ /• 

Jm? fc2 + S Jm^ ^ Jrt 



s A:2 + s 



+~ , m2 1 

ds 

,2 s s 



+°° , m2 1 /•+°° , m2 1 

cts — 7^^— h / ds 

s k-' + s J ^2 s s 



= J_ r°° as 

167r2 



_A;2 + s s 

We conclude that the spectral density is given by 



which is indeed positive for s> m^. 



Appendix C: Derivations of the spectral density for the simplest correlation functions in 

Gribov-Zwanziger theory 

1. Using the Feynman parametrization for d = 2 



We wish to bring 



31 



into a spectral form, with F{k'^) defined in eq. ( 120 ). We first derive a Feynman parametrization of ( 120 1. Proceeding 
in the usual way one finds 



d-^q N{q,k,x) 

+ A2)2 ' 



where we used the substitution q = p — kx, and whereby 

= x{l - x)k^ - {2x ~ 1)^7^ 



(C2) 



(C3) 



We shall temporarily work in units 27^ = 1. We still have to identify the numerator N{q, k, x). Keeping in mind that 
terms odd in q^ will vanish upon integration, and that we may replace g^t?!/ — > 9^%^ within the g-integral, we are 
brought to 



F(fc^) = {d-l) dx 



1 '• A x^{l-xfk^+l[l-{d + 2)x{l-x)]k^q^+q'' 



after a bit of algebra. We then recall that 

d'^q 1 



r d^q 
J (27rF 



(92 + A2)2 
1 T{72 - d/2) 



(g2 + A)" (47r)^/2 V{n) 



(A2)'*/2- 



from which it follows that 



and consequently also 



d'^q 



i <i r(n-rf/2- l) ^^2-)rf/2-n+l 



{2'kY (g2 + A)" (47r)''/2 2 T{n) 



d'^q 



t ^ 1 d{d + 2)r{n-d/2-2) , 

(27r)'' ((j2 + A)" (47r)'^/2 4 r(n) ^ ^ 



To obtain a finite result, we prefer to look at 



(9fc2)2 - ^ 



dx 



1 

A2 



(I2a;'^ - 24x3 _^ ^4^2 _ 22:) — + (8.t^ - 24a;^ + 25a:'' - lOa;^ + x^^ 



(A 



2^2 



+2x^(1 - xy 



(A2)3_ 

where we set d = 2, as we are mainly interested in this case now. 
We consequently find 



~(dk 



2^2 



4^ 7o 



-12x2 + 12x-2 ,8x2-8x + l , x(l - x) 
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where we reintroduced s — 2x{i-x) ' ti^nce x = ''"''"''t/"'" "^; which gives rise to 



{dk^f 



1 

An 



+00 



2s 

ds 



2(1 + s2 + Vl + s2) 

^(4 + s2_4yrT^)-^ 



K3 + .s2-3v/l + .s2)— ^ 
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We first rewrite everything in terms of k^ — is as follows 

ds 
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(C6) 



(C7) 



(C8) 



(C9) 
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1 
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(Cll) 



Using two consecutive partial integrations, we can show that 



2(1 + s2 + yrr^) 
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Similarly, partial integration leads to 
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(C12) 
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Hence, we can rewrite (Cll I as 



92f (fc2) _ ^ 
(9fc2)2 ^ 4^ 



+00 



ds 
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(C13) 



(C14) 



after simplification. We observe that there are no poles in the upper half s-plane for k"^ > 0, so we can fold our 
contour around the cut for s G [ioo, i]. With s = ir, we can write 



92F(fc2) _ ^ 
47r 
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idr 
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We can now return to the original function F{k'^). A first integration from to k^ gives 
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so that we get 
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Next, we shall also bring G(fc2), give in eq. (1251, into a spectral form. Invoking the Feynman trick this time yields 

r-l 



G(fc^) 



dx 



{2TrY (g2 + A2)2 ' 



with 



M{k,q) = k^q^ 1 



As usual, we shall work in units 27 — 1. By using (C6), we obtain the finite function 
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4J\f I dad/?cos(f (a-/3)) I {p''{p - k)^ + {d ~ 2){p^ - p kf) e-^P'-^^^-'^^" 



(27r 



= 32A/' / da dp cos (7^ (a - 
Jo 

The Gaussian integrals can be evaluated by standard methods: 

/•oo 

(0S(fc)0S(-fc)> = W / dad/? cos(7'(a-/3)) (d- 1) 



34 
(C31) 

(C32) 



(0S(fc)0S(-fc)>=32AA / dadpcos{f{a-p)){d-l)Ua + P)-'"'-' 



(C33) 
(C34) 



For the remaining calculation it is useful to obtain a general expression for such integrals as appearing above. They 
can be derived along the lines of Section IIIB only that we consider now arbitrary powers of a, /3 and {a + j3). 



We define the expression J as 







and insert unity by means of eq. (A2|. Then we rescale the variables a and /3 by A, 

»i 
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evaluate the integrals over A and /?, and transform a to the variable u = (2a — 1)/(2q!(1 — a)): 

11 1 
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This integral is evaluated in the complex plane, where it possesses a branch cut from u = z to m = ioo. Setting u — iy 
the integral J is then given by 



f°° 1 
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(C38) 



Now the physical cut from fc^ — —272 to infinity is manifest, but we still have to determine the spectral function. For 
this we rewrite the variable y into an angular variable x- — cosx, < x ^ ""/S, 



J(a, b, c, d) = (47r)-'^/2 2^+''T{2 + a + b + c)i dx 
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.77/2 C 



(C39) 



and employ eq. (AlO) to obtain 

J{a,b,c,d) = (47r)-'^/2 22+c 



dr- 



T{-l-a-b-c) J2^2 T + k^ 



(C40) 



where the order of integration was interchanged and ip is defined by cos (p = 26'^ /t. Note that the external momentum 
k'^ in a coeflicient of J will turn into — r. 

Eq. ( C40 ) is valid for —2 — a — 6 — c 7^ —1, —2, .... In our example there is indeed one case, where —2 — a — b — c = — 1 
is realized. Here we use instead the intermediary expression eq. (C39l of the deriv ation above . It leads direc tly to 
the spectral representation after replacing cos(x) by Finally plugging eqs. (C40l and (C39l into eqs. (C33l 

and ( C34 1 yields in two dimensions 



-,(1) 



dr 



12^2 T + k^ 271 ./t^ - 474 
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These results are consistent with (1231 and (1281. In four dimensions, we find 



1 -474(274 + t2) 



32tt^t 

3/2 



(C41) 
(C42) 

(C43) 
(C44) 



[1] V. Mathieu, N. Kochelev and V. Vento, Int. J. Mod. Phys. E 18 (2009) 1. 
[2] V. N. Gribov, Nucl. Phys. B 139 1 (1978) 
[3] D. Zwanziger, Nucl. Phys. B 323 (1989) 513. 
[4] I. M. Singer, Commun. Math. Phys. 60 (1978) 7. 

[5] G. Deir Antonio and D. Zwanziger, Commun. Math. Phys. 138 (1991) 291. 

[6] D. Zwanziger, Nucl. Phys. B 412, 657 (1994). 

[7] D. Zwanziger, Nucl. Phys. B 399 (1993) 477. 

[8] N. Maggiore and M. Schaden, Phys. Rev. D 50 (1994) 6616. 

[9] D. Dudal, R. F. Sobreiro, S. P. Sorella and H. Verschelde, Phys. Rev. D 72 (2005) 014016. 
[10] J. A. Gracey Phys. Lett. B 632 (2006) 282. 
[11] J. A. Gracey JHEP 0605 (2006) 052. 

112] D. Dudal, J. A. Gracey S. P. Sorella, N. Vandersickel and H. Verschelde, Phys. Rev. D 78 (2008) 065047. 
[13| L. Baulieu and S. P. Sorella, Phys. Lett. B 671 (2009) 481. 
[14| S. P. Sorella, Phys. Rev. D 80 (2009) 025013. 

[15| D. Dudal, S. P. Sorella and N. Vandersickel, arXiv: 1001. 3103 [hep-th[. 

[16| R. Alkofer, W. Detmold, C. S. Fischer and P. Maris, Phys. Rev. D 70 (2004) 014014. 

[17| K. Osterwalder and R. Schrader, Commun. Math. Phys. 31 (1973) 83. 

[18| K. Osterwalder and R. Schrader, Commun. Math. Phys. 42 (1975) 281. 

[19| L. von Smekal, R. Alkofer and A. Hauck, Phys. Rev. Lett. 79 (1997) 3591. 

[20| L. von Smekal, A. Hauck and R. Alkofer, Annals Phys. 267 (1998) 1 [Erratum-ibid. 269 (1998) 182]. 

121] D. Zwanziger, Phys. Rev. D 65 (2002) 094039. 

122] C. Lerche and L. von Smekal, Phys. Rev. D 65 (2002) 125006. 

[23] J. M. Pawlowski, D. F. Litim, S. Nedelko and L. von Smekal, Phys. Rev. Lett. 93 (2004) 152002. 
[24[ M. Q. Huber, R. Alkofer and S. P. Sorella, Phys. Rev. D 81 (2010) 065003. 



36 



[25] J. C. Collins, Renormalization. An Introduction To Renormalization, The Renormalization Group, And The Operator 

Product Expansion, Cambridge Univ. Pr. (1984). 
[26] D. Dudal, S. P. SorcUa, N. Vandcrsickol and H. Vcrschcldc, JHEP 0908 (2009) 110. 
127] G. 't Hooft and M. J. G. Veltman, Scalar One Loop Integrals, Nucl. Phys. B 153 (1979) 365. 
128] C. Itzykson and J. B. Zuber, Quantum Field Theory, McGraw-Hill (1980). 
129] D. Dudal, S. P. Sorella, N. Vandersickel and H. Verschelde, Phys. Rev. D 77 (2008) 071501. 
[30] D. Dudal et al., work in progress. 

[31] E. Seilcr, "Gauge Theories As A Problem Of Constructive Quantum Field Theory And Statistical Mechanics," Lect. Notes 

Phys. 159 (1982) 1, Berlin, Heidelberg, New York; Springer (1982). 
[32] K. Osterwalder and E. Seller, Armals Phys. 110 (1978) 440. 
[33] M. Luscher, Commun. Math. Phys. 54 (1977) 283. 
[34] D. Zwanziger, Nucl. Phys. B 378, 525 (1992). 

[35] N. N. Bogolyubov, A. A. Logunov and I. T. Todorov, "Introduction To Axiomatic Quantum Field Theory", Reading (1975). 
[36| U. Habel, R. Konning, H. G. Reusch, M. Stingl and S. Wigard, Z. Phys. A 336, 423 (1990). 
[37| U. Habel, R. Korming, H. G. Reusch, M. Stingl and S. Wigard, Z. Phys. A 336 (1990) 435. 
[38| M. Stingl, Z. Phys. A 353, 423 (1996). 

[39| L. Driesen, J. Fromm, J. Kuhrs and M. Stingl, Eur. Phys. J. A 4, 381 (1999). 

[40| L. Driesen and M. Stingl, Eur. Phys. J. A 4, 401 (1999). 

[41| M. Stingl, Phys. Rev. D 34, 3863 (1986) [Erratum-ibid. D 36, 651 (1987)]. 

[42| G. V. Dunne and I. G. HalUday Phys. Lett. B 193 (1987) 247. 

[43[ R. Ricotta, in J.J. Giambiagi Festschrift, ed. H. Falomir (1990) 350. 

[44| I. G. Halliday and R. M. Ricotta, Phys. Lett. B 193 (1987) 241. 

[45] G. V. Dunne and L G. HaUiday Nucl. Phys. B 308 (1988) 589. 

[46] F. A. Berends, M. Buza, M. Bohm and R. Scharf, Z. Phys. C 63 (1994) 227. 

[47| S. Bauberger, F. A. Berends, M. Bohm and M. Buza, Nucl. Phys. B 434 (1995) 383. 

[48[ C. Anastasiou, E. W. N. Glover and C. Oleari, Nucl. Phys. B 572 (2000) 307. 



